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Abstract 

The Weyl-Wigner formulation of quantum confined systems poses several interesting problems. The 
energy stargenvalue equation, as well as the dynamical equation does not display the expected solutions. 
In this paper we review some previous results in the subject and add some new contributions. We 
reformulate the confined energy eigenvalue equation by adding to the Hamiltonian a new (distributional) 
boundary potential. The new Hamiltonian is proved to be globally defined and self-adjoint. Moreover, 
it yields the correct Weyl-Wigner formulation of the confined system. 
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1 Introduction 

Deformation quantization is a powerful procedure to derive the quantum mechanical formulation of a wide 
range of physical systems. These include systems defined in non-flat Poisson or sympletic manifolds [1]. 
In the simplest case of flat phase spaces the deformation approach yields the Weyl-Wigner formulation of 
quantum mechanics [2]- [9]. This is an autonomous formulation of quantum mechanics that has been an 
active research topic in both physics and mathematics. 

Quite unexpectedly the deformation approach displays several problems in the case of confined systems 
[10]- [13]. The dynamical and the eigenvalue equations do not yield the correct solutions. Given the status of 
the deformation method this is a serious drawback. Moreover the importance of confined systems in several 
domains of mathematical physics [14]-[16] (e.g. the theory self-adjoint extensions of symmetric operators), 
condensed matter physics (e.g. the quantum description of particles moving on surfaces with obstacles or 
impurities) and even in string theory and other modern approaches to quantum gravity [17] (where the 
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classical theory displays a non trivial topological structure) justifies a more thorough investigation of the 
quantum phase space formulation of these systems. 

In this paper we study the simplest version of this problem: the Weyl-Wigner formulation of the en- 
ergy eigenvalue equation for a 1-dimensional particle confined to the negative semi-axis and subject to 
Dirichlet boundary conditions. We will identify the source of the problem at the level of the standard oper- 
ator formulation of quantum mechanics by showing that the confined energy eigenfunctions are not global 
(i.e. distributional) solutions of the textbook eigenvalue equation. We will propose a new, globally valid 
eigenvalue equation which is written in terms of a new Hamiltonian containing a confining distributional 
potential. The translation of this equation to the Weyl-Wigner formulation yields a new stargenvalue equa- 
tion displaying the correct solutions. This problem has been studied before [10]- [13]. Here we reformulate 
it on simpler grounds and add some new contributions. In particular we show that the new Hamiltonian 
is self-adjoint. From the point of view of the theory of self-adjoint extensions of symmetric operators this 
is also a non-trivial result [14]- [16]. 

This paper focuses on the stargenvalue equation but a similar approach is valid for the dynamical 

2 

equation. The method presented here is extendable to generic Hamiltonians of the form H = ^ + V(q) 
and to general boundary conditions. 



2 The Weyl-Wigner formulation of confined systems 

Let us consider the Gelf'and triple Q C 7~i = L2(IR) C Q' where Q is the space of infinitely smooth functions 
t(x) that, together with all its derivatives, decay to zero faster than any power of 1/x as \x\ — > oo, Q' is its 
dual and L2OR) is the space of complex valued square integrable functions on IR. Generalized operators are 
linear maps from Q to Q' . The set of these operators is denoted by T and (in a suitable sense) it contains 
the sets of Hilbert-Schmidt operators over TC as well as the enveloping algebra of the Heisenberg-Weyl Lie 
algebra. The Weyl-Wigner map is a one to one invertible transformation from T onto the set of generalized 
functions on T*M ~ IR 2 (we make h = 1) [5]: 

A — > A(x,p) = W[A}= J dye-' ipy < x + y/2\A\x - y/2 > . (1) 

This map provides the translation of all the important structures of operator quantum mechanics into 
the quantum phase space language: operators are mapped to phase space functions (denoted by symbols) 
A(x,p) = W(A), and the density matrix to the celebrated Wigner quasi-distribution 

f w (x,p) = ^W(\^j X VI) = f dye'*™ <x + y/2\i> >< ^\x - y/2 > . (2) 

The operator product and bracket yield the algebraic structures of the non-commutative phase space. We 
have 

'i ( d d d d 



A(x,p) * B(x,p) = exp 



A(x,p)B(x',p'] 



(3) 

(x',p')={x,p) 



2 \ dx dp 1 dp dx' , 

for the ^-product and [A, B\ M = —i(A * B — B * A) for the Moyal bracket. The two key equations of 
the Weyl-Wigner formulation are the Moyal dynamical equation and the *-genvalue equation. The former 
reads 

A * F aa = F aa * A = aF aa . (4) 

It amounts to the Weyl transform of the eigenvalue equation of the standard operator formulation and is 
of equal importance in the context of quantum phase space methods. Its solutions (the *-genfunctions) 
are ^.-deformations of the Dirac measure, being central to the calculation of probability distributions. 
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Let us now consider a one dimensional quantum particle described by the Hamiltonian (H = 2m = 1) 



H =p 2 



8P_ 

dx 2 ' 



(5) 



which is confined to the negative semi-axis x < and subject to Dirichlet boundary conditions at x = 0. 
The energy eigenstates of the system are 1 



^e(x) = 9(-x)Mx) = 0(-x) 
The Weyl-Wigner transform of which are 

Fee(x,p) = 6{—x 



E>0. 



2sin ^ (p+ ^] + 2s ^-^>] + | cos( 2xv^) sin(2xp) 



(6) 



(7) 



Contrary to our expectations these phase space functions do not satisfy the energy stargenvalue equations 
(not even for x < 0) 

H*F EE ^EF EE and F EE *H + EF EE . (8) 

The questions we want to answer are the following: (i) Why aren't the stargenfunctions solutions of 
the stargenvalue equation? (ii) Is it possible to modify the Hamiltonian so that it generates the proper 
solutions? (iii) Is the new Hamiltonian self-adjoint? 

To look for some guidance let us go back to the operator formulation of quantum mechanics. Substi- 
tuting the eigenfunctions (6) in the eigenvalue equation we easily find that they are not global eigenstates 
of (5). In fact 



HoM*) = [8{-x)4>e{x)} = E^e(x) + 25{x)c/> E {x) + 5'(x)<p E (x). 



(9) 



The new distributional term on the rhs is not zero for Dirichlet boundary conditions. Hence, the energy 
eigenvalue equation is also not globally valid for confined systems. Since the Weyl-Wigner map applies to 
the set of generalized operators it comes as no surprise that eqs.(8) are not valid either. 



3 Global formulation of the quantum confined particle 

We now propose an alternative formulation of the confined energy eigenvalue equation which: (i) Is valid 
in distributional sense, i.e. the confined eigenfunctions are generalized solutions of the new eigenvalue 
equation, (ii) is formulated in terms of a Hamiltonian displaying a new confining boundary potential, and 
such that (iii) the new Hamiltonian is self-adjoint. 

The first step is the introduction of the following definitions: 
Definition 1: Special distributions of order n 

A Schwartz distribution 2 F is a special distribution of order n iff for each xq £ HI there is a left and a 
right neighborhood of xq, =\xq — o"_,xo[ andVL + =]xo,XQ + a + [, (o~-,cr + > 0), and two C n (IR) -functions 
/_ and f + such that, in the sense of distributions, "F = /_ on Sl„ " and "F = f + on Q + ". The set of 
points where F is not a C n (IR) -function is denoted by sp supp n F and the set of special distributions of 
order n by S n . In particular we have sp supp 00 ^ =sing supp F [18]. 

This definition provides a hierarchy of spaces 5°° C ... C S 1 C S°. Trivial examples of elements of <S°° 
are the Dirac delta 5 and all its derivatives and of 5° are all limited functions which are continuous except 

4>e(x) are the free eigenfunctions and 6(x) is the Heaviside step function. 
2 The space V of Schwartz distributions [18, 19] is the dual of the space V of infinitely smooth functions of compact support. 



3 



on a finite subset of IR. 

Definition 2: The regularized delta operators 5±\x — a), ogE 
5^ (x — a) are the linear operators 

5 { ±\x-a) : S n ^S°°; 

F — ► 5 { ±\x - a) [F] = lim 5 {n \x - a) ■ F(x ± e). (10) 

where the product ■ appearing on the rhs is the one proposed by L. Hormander in Ref.[19j. 

Notice that for e g]0, a[ and a sufficiently small we have sing supp 5^ n \x — a)f\ sp supp n F(x±e) = {a}n 
sp supp n F(x ± e) = 0. Hence the product of distributions on the rhs of eq.(10) is well defined in the sense 
of the product of distributions with non-intersecting special supports proposed in [19]. This product yields 
o~(n)( x — a)f±(x ± e), where f± is the C™-function associated to F in the right (left) neighborhood of a (cf. 
Definition 1). Consequently, the limit of (10) exists in 5°°. 

Let us illustrate the action of these operators with two simple examples 

5 ( +\x) [e^ix)] =5 mfi 5 ( - n) (x) and (x) (x)] = 0. (11) 

Using the operators 5±^ we can re-write the energy eigenvalue equation (9) as 

~^2^e{x) = - [0(-x)Mx)] = EiPe(x) + 26-(x)il/ E (x) + 5'_{x)^ E {x). (12) 

Notice that for ip(x < 0) = <p{x < 0) and 4>(x) G C°°(]R) we have, in the distributional sense, 25- (x)ip'(x) + 
S'_(x)i/j(x) = 25(x)4>' (x)+5' (x)4>(x). Imposing the Dirichlet boundary conditions ?Ae(0) = 0, eq. (12) reduces 
to 

-^e(x) = E^ E (x) - 5'_{x)i> E {x). (13) 
Hence, the new (confining) Hamiltonian is 

H-f + tM.tk + tM. (14) 

2m 2m 2m 

We can solve the new eigenvalue equation (13), both analytically and numerically (with a suitable smooth 
regularization of the Dirac delta) [10], for Dirichlet boundary conditions and show that all its solutions are 
of the form (6). 

To put this formulation on firm mathematical grounds, we now specify the domain of the Hamiltonian 
and prove that it is self-adjoint. Let us assume the Hilbert space to be 7i c = {ip € L2(JR) : ip(x) = 
9(—x)(/)(x)} and consider the action of H. Its maximal domain (i.e. the largest domain where H is well 
defined) is 3 

T> rnax = {4;eH c :4> = 9(-x)cf>(x), 4>{x) G AC 2 (M), 0(0) = 0}. (15) 
To prove this let us consider the action of H on a generic if>(x) = 9(—x)(f)(x). We have 

Hif,(x) = - [6{-x)4>"{x) - 25{x)(t>'{x) - 5'{x)4>{x)] + 5'_(x)ip(x) 

= -6(-x)<f>"(x) + [25(x)<f>'{x) + 25'(x)<t>(x)] , (16) 

and so 

Hij)(x) G H c => 5{x)<p'{x) + 5'{x)4>{x) = =^ ^(0) = 0. (17) 



AC (IR) denotes the set of functions defined in IR with absolutely continuous first derivative. 
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Hence H acts as 



H:V n 



H c ; H [9(-x)(t)(x)] = -0(-x)<f>"(x). 



(18) 



Notice that V max is dense in 7i c and that H is closed in V max . To study the self-adjointness of H let us 
introduce the sesquilinear form 



and we easily find that 



+00 



MM) = e (o)v^(o) -^(O)^'(O) = 0, V£,V G P„ 



(19) 



(20) 



and so H is symmetric in its maximal domain. Therefore H is also self-adjoint. 

To check this result we may use von Neumann's theorem on deficiency indices. For x < the solutions 

of 

(H ± ik 2 )^j = 0, k>0, (21) 



which belong to H c , are Y>±(x) = C± exp ^^ L x and none (except for the trivial solution) satisfies the 

Dirichlet boundary conditions at x = 0. Hence, the deficiency indices of are (0, 0) meaning that H is 
indeed self-adjoint. 



4 Wigner formulation of confined systems 

Using the Hamiltonian (14) we may now write the new stargenvalue equation 



2m 2m y ' 



★ Fee = EFee and Fee * 



f- + ±fi (x) 
2m 2m y ' 



= EFee, 



where the distributional term is evaluated according to 



5'_(x)*F E e(x,p) = lim [S'(x)*F E e(x - e,p)] , 



and impose the boundary conditions: 

V^(0) = ^MO)Ve(O) = /" d P F EE (0,p) = 0. 



(22) 



(23) 



(24) 



Equation (22) with the boundary conditions (24) displays the unique solutions (7). This was proven in 
[10]. 
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